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^^ ' 1. Introduction 

^ I In this paper, we revisit some results by Nieto and Truax about generating functions for arbitrary order coherent 

and squeezed states. These results were obtained using the exponential of the Laplacian operator; they were 
later extended by Dattoli et al. [8J using more elaborated operational identities. In this paper, we show that 
OQ ■ the operational approach can be replaced by a simple probabilistic approach, in the sense that the exponential of 
derivatives operators can be replaced by equivalent expectation operators. This approach brings new insight about 

r^ . the links between operational and probabilistic calculus. 

O ' In the first part, we show that the exponential of the derivation operator of arbitrary integer order can be replaced 

^ . by an expectation with repect to a carefuly chosen random variable. 

In the second part, we apply this result to the Gould-Hopper polynomials, which include as special cases the 
Kampe de Feriet and the Hermite polynomials; this allows us to recover and interpret easily some properties of 
these polynomials. 

C/3 . The last part is dedicated to the application of these probabilistic representations to the computation of coherent 

"*-i ' and squeezed states. 

s' 

'^ ■ 2. The operator exp 

O . 2.1. Introduction. In j2j, Nieto and Truax consider the operators 



a 



(N 

in 

o 
o 



d \1 



dx 



(SI ■ ^j = exp 

> 

vQ ' where c is a constant and j an integer. They remark that, for any well-behaved function/, /i acts as the translation 
operator 



in ; 

^P I Moreover, Z2 being a Gaussian random variable with variance 2, I2 acts as the Gauss- Weierstrass transform 

hf (x) = EzJ (x + CZ2) . 



Since Iif can be also written as Ezif (x + Zi) where Zi is the deterministic variable equal to 1, it is tempting to 
K% ' wonder if the expression 

. .C^; I, ^ EzJ {x + cZ,) 

holds for values of j > 3 and to study what random variable Zj possibly comes out in this formula. In [2-, the 
following general expression is proposed 



(2.1) I, fix) 



1 f+°° da;ie-^i /■+°° dxafi-^^ f+°^ dx.^iC''-^ 



X 



{271)' ^'^jh x\ J a ^1 Jo ^^J 



and the authors add: "This result, although in closed form, is usually too complicated to yield results in terms of 
elementary functions". The next section is devoted to providing a simple expression for the operator Ij. 
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2.2. a complex valued stable random variable. We first introduce the following random variables Wj and Zj 
where in this whole section, j is an integer such that j > 2. 

Definition 1. For j G N, j > 2, define the root of the unity 

/ 2n 
Wj = exp I i^ 

and the random variable Wj as 

(2.2) Pr {VF = wj} = i, < ; < j - 1. 
Thus Wj equals equiprobably all values of the j— roots of the unity. The random variable Zj is then defined as 

Zj^W.Sl'^ 

j 

where Si is a stable random variable [6J with characteristic parameter - , independent of Wj . 

In the following, when possible, we'll discard the indices for notational simplicity and write simply 

(2.3) Z^WS'T. 
We now describe some properties of the random variable Z. 

Proposition 2. The random variable Z has moments 



(2.4) EZ'' = 

p 

Proof. The k ~ th order moment of Z is 



M! tfk^pj,peN. 



with 



EW" 



Moreover, with k ~ pj, 



EZ^ = 


^EW'' 


Si 








1 ^"^ 

■' 1=0 
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ifk^ 
if fc = 


PJ 
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(PJ)! 
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(see [3]) so that the result follows. D 

As a consequence, the generating function ipz (u) = E exp {uZ) of Z can be evaluated as follows. 
Proposition 3. The generating function of Z is 

'■Pz (u) = exp (u-') , u> 0. 
Proof. A straightforward computation gives 



k=0 ■ p=0 ^-^ ■ ^' 



1\ 



n 

From this property, we deduce that the random variable Z exhibits the following stability property. 

Proposition 4. // Z\ and Zi are two independent random variables distributed as in i2.3\) with parameter j and 
if fli and 02 are two real positive numbers then the random variable 



aiZi + 02^2 ^ (a{+ 0^2)' Z 

where Z is again distributed as in 112. 3\) with same parameter j and ^ denotes equality in distribution. 

We note that Z is a complex valued random variable; for real valued random variables, the stability property 
can hold only for values of the parameter j such that < j < 2 (see ^ p. 170]). 
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2.3. application to the operator exp 

follows. 



d \1 



j^) . We can now give a simplified version of the formula p.l|) as 



acts on a function f as 



Theorem 5. The operator Ij — exp 

(2.5) I,fix)=EzJ{x + cZ) 

where Zj is defined as in 112. 3\} . for any function f such that the right-hand side of this equality exists. 

Proof. By definition, with Z defined as in 



Ijf (x) = exp 



dx 



f (x) = Ez^ exp ( cZ— ]f{x)^ Ez^ f [x + cZ) . 



n 



We note that the link between equalities (|2.ip and (|2.5p can be identified using Williams' formula |3l end of part 
2]: formula (|2.ip can indeed be rewritten as 



Ijl {^) = EwEy^.y 

3 

where Yi , Y2 , . . . , Yj-i is the product of j — 1 independent Gamma random variables with respective shape param- 



eters 



1 2 



. , '^-J-. But by Williams' formula, 



\j J 1 
where Si is a stable random variable with parameter -, so that 

I,f{x)=Ew^Sif{x + cWS7^ 



which is exactly (|2.5[) . 

As an example, we consider the case j = 2 for which W is Bernoulli distributed with parameter i and since Si 
is Levy distributed, according to f¥5 p. 456], 

^ ^ V2|7V| 



where N is Gaussian with variance aj.^ = 1 so that Z2 — W2S1 ^ is itself Gaussian with variance (t| = 2 and we 

2 
recover the Gauss- Weierstrass transform operator 

hf{x)^Ef{x + cZ2). 

3. Properties of the Gould Hopper polynomials 
As an application of the former result, we consider the particular case / [x) = x". The polynomials 



(3.1) 



9n {x,h) = exp 



dx 



x", {m,n) e 



are known as the Gould-Hopper polynomials, as introduced in [T] p. 58]. The particular case m ~ 2 corresponds to 
the Kampe de Feriet polynomials, and the polynomials g^ {2x, —1) coincide with the classical Hermite polynomials. 
From the results of section [2] and with h — c"^, we deduce 

(3.2) 



gl^{x,h) = Ez[^x + h-Z 

where Z is distributed as in (|2.3I) with j = m\^. We note that this representation is different from the one given by 
[5j Thm.l], which seems of little use since it involves integer moments of stable random variables. 

We show now that the representation p.2p allows to recover and extend easily some well-known properties of 
the GH polynomials. 



we replace here j by m to stick to the notation introduced by Gould 



A PROBABILISTIC APPROACH TO SOME RESULTS BY NIETO AND TRUAX 



(1) the generating function of the GH poynomial is [TJ (6.3)] 



n=0 n=0 



= Ez cxp (tix + h^zYj = exp {tx + ht"") ; 

(2) the derivative of a GH polynomial [U (6.4)] is easily computed as 

—5™ (x, h) = nEz [x + h~Z\ = ng™_^ (a;, h) . 

(3) the GH polynomials satisfy the following addition theorem 

(r r \ r 

■i=l 1=1 / riiH hnr=n/c=l 

the proof of which is straightforward using the representation (|3.2p . This result generalizes the case r = 2 
obtained by Gould and Hopper under the form fT. (6.19)] 

E it) 5" (^' '^^ 5-^ (y' '^^ = 5" (^ + y, 2/1) . 

fc=0 ^ ^ 

Moreover, the classical addition theorem for Hermite polynomials [7, (40) p. 196] 

is easily deduced from (|3.3p using the link between Hermite polynomials and Gould polynomials mentioned 
above. 

(4) Another result from [1, (6.9)], namely 

exp [h-^ j (a;" exp {tx)) = — (exp {tx + ht^)) 



can be recovered easily using (|2.5p since 

exp I /i-— I (x"exp(ta::)) = Ez (x + hi Zj exp (t (x + h^ Z 



= rf^^^^^^P (^ (^ + '^^^j j ^ d^ («^P (^=^ + ^*')) ■ 

4. Application to Coherent and Squeezed States 

4.1. Multisection of series. In calculations of higher-order and squeezed states quantities, the following gener- 
ating functions of Hermite polynomials appear - we follow here the notation by Nieto and Truax: 



(4.1) S{j,Kz) ^ Y. 



jn+k 



^-^ {jn + k)\ 

n=0 ^-^ ' 

+ °° ^jn+k TT ( ^\ 



n=0 

and 

+°° ^jn+k 

(4.3) g{j,k,x,y,z) = Y 



{jn + ky. 
z^"+'=7Ij„+fe {x) H,r.+k {y) 



(in + k)\ 
n=o ^-' ' 

where _ff„ {x) is the Hermite polynomial of degree n, and j and k are integers. 

These quantities are all computed in [2 using an operational calculus approach, but we show here that the 
probabilistic approach introduced in Section [5] gives a new insight on these formulas and allows to extend them to 
other classes of polynomials. 

We introduce the following theorem. 
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Theorem 6. // 

+ 00 

n=0 

is the generating function of the sequence of functions {/„ (x)} then]j 

+00 

where W is distributed as in i2.S\). 



Proof. We compute 

+00 

Ew,W-''cp{tW,x) = ^Svy^.W^"-'=i"/„(x) 

and use the property p.4p to deduce the result. D 

We note that this theorem is the probabihstic formulation of the technique of "multisection of series" described 
by Riordan [SI section 4.3] as a "process of ancient vintage". This process is also characterized in [8, (19)] as a 
consequence of the sieving principle. 

4.2. Applications. As an application of this theorem, we recover easily the results in [2]: 

(1) choosing /„ (x) — -^ Vx gM., Vn G N we deduce 

(j> (t, x) = cxp (t) 

and 

S ij, k, z) = Ew, W-'' exp (zW) 

(2) choosing /„ (x) — ", -^ , we deduce the well-known generating function of the Hermite polynomials 

(j){t,x) =exp(2te-i^) 

so that 

G (j, k, X, z) = Ew, W-^ exp {2Wxz - z^W^) 

(3) the more general case 

+00 pn+k 
n=0 ^•' ' 

can also be easily derived: choosing /„ (x) = "+'" 1 we deduce the generating function ^llj (1) p. 197] 

(j) {t, x) = exp (2ix - t'^) H„i (x - t) 
so that 

+°° fjn+k 

J2 (■„ , M! -^J»+fe+'» (^) = ^M'. ^"^ exp (2Wa; - t^W^) i/„ (x - W) 

n=0 ^■^ '■ 

which coincides with [8, (21)]. 

(4) remarking that the variable x in theorem ^ may be multidimensional, choosing x = (xi,X2) and 

Hn{xi)Hn (X2) 

/ a; = j 

we deduce, by the Mehler formula [II (22) p. 194], 

1 / t'^x^ + t^'iP' — txy 
Ht,x,,x,) ^ -^==exp (^-4 ^-^^^ 

so that 



1 / t^W^x^ +t'^W^y^ -tWxy\ 



q(j,k,x,y,z) — EwW ' , exp —4 ; . „„^„ 

We remark that Theorem [6] is not restricted to Hermite polynomials; let us give another few examples: 



^we recall the notation Ew^f {W) = i ZlLo f (^' 
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(1) for the Gould Hopper polynomials g™ (x, h), choosing /„ (x) = ^" , ' as defined in p.ip . we deduce 

+°° pn+k 

Y. T—rV\\9rn+k (x, h) = Ew, W-" exp (txW + ht^W"^) 



^ (jn + ky.- 



r!=0 



^"^1 ^^r^ I xt 



(2) for the Laguerre polynomials, with /„ (x) — L^ (x) we have (f>{t,x) = (1 ^ ^) cxp ( -^ J so that 

+2? j-jn+k / ^f]^r 



4.3. Further sums. In their study, Nieto and Truax consider also the sums 



+ 00 



zJ"+'=i/,„+fe [x) 



(4.4) X(j,fc,p,g,x,z) = ^ 

for which they explicit some specific cases. We give here the general result only in the cases j — p and j — 2p by 
looking first at the simple sum 



+°° ^jn+k 



(4-5) k{j,k,p,q,z) = V , ... 

A straightforward computation shows that this sum is related to the sum S (j, k, z) in (14. ip as 
(4.6) k{j,k,p,q,z) = z''^^pS[p,q,zij . 

As a generalization of this result, we obtain the following 
Theorem. If the sequence of functions {/i„ (x)} can be expressed as a sequence of moments 

K {x) = Eu^'' {x, U) 
for some function Lp, then 



(4-7) E J„_^g), V^+fc {x) = i?w,,t/ [{zv {x, U)f-'i Ty-«exp ((z^ (x, t/))^ M^) } . 

Proof. The proof is straightforward replacing z by z(y9 (x, U) in (|4.6p and taking expectation over U. D 

As a consequence, using the moment representation of the Hermite polynomials 

Hn {x) = 2"^Ar (a; + iN^- 
where A'^ is a Gaussian random variable with variance cr^ = i we deduce 

(4.8) K (j, k,p, q, x, z) = {2zf^^i Ew^.nW''' {x + iNf~^i exp ({2z {x + iN))i w] 

The expectation with respect to the Gaussian variable N is difficult to obtain for an arbitrary value of the ratio - 
so that in the following, we give explicit expressions only in the cases j = p and j — 2p. 

4.3.1. the case j ~ p. In this case we obtain the following result. 

Theorem 7. // fc — q G N, the function K {j, k, j, q, x, z) reads 

K{j,k,j,q,x,z) = z^~'iEw,W-'ieyi^{-zW{zW -2x))Hk-q{x-zW) 

Proof. The proof is immediate remarking that in (j4.8p the Gaussian expectation 

ENix + iNf^''eyi^{2zW{x + iN)) = (21^)*"'' -— -^ exp (2x2^ £;Ar exp (i2zM^Ar) 

dz'^~''i 

involves the Gaussian characteristic function i^Arexp (i2zM^A^) ~ exp (— z^PI^^) so that it is equal to 

(214^)'-'= exp (x2) ^ exp (- {zW - xf) 
which, using Rodriguez formula for the Hermite polynomials, can be expressed as 

2'^-^ exp {x^) exp (~ {zW - xf\ {-if^'' Hu-q {zW ~ x) 
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from which we deduce 

K (j, fc, j, q, X, z) = z'^-'iEw, W^'i exp {-zW {zW - 2x)) Hk-q {x - zW) 

4.3.2. case j — 2p. This case can also be solved as follows. 
Theorem 8. If k — 2q E N, the function K {2p, k,p, q, x, z) reads 

exp (jj^) 
Ki2p,k,p,q,x,z) = z^-^'^Ew.W-'i \ \_l^, Hk-2q 



n 



Proof. With j = 2p, replacing x by ix, we need to compute the Gaussian expectation 

En {x + Nf~'^'^ exp i-Az'^W {x + Nf\ = —= j (x + Nf~'^'' exp i-'iz^W {x + Nf\ exp {-N"^) dN 



which, after change of variable M — x + N reads 

"+00 



1 /"^ 

—= exp {~x^) / A'f '=^-2? g^p |-_^^^2 (^^ ^ /iz^w)) exp (2xAf ) dM. 

V"" J-00 

This classical Gaussian integral is equal to [10, 2.3.15.10] 

- exp [-~x ) T-^+T Hk-2q 

.2/ (l + 4z2t^)^^ \\/l + Az^W 

and we deduce, replacing x by (— ix), 

cxD f "^^'y xA 

if (2p, A:,p, q, X, z) = z'^-^'^Ew.W-^ ^- ,r^Hk-2q 



k-2g+l --/i-iiy 1 I- ; — - — yYjy 

(l + 4z2W^) — ^ — KVI + Az^W 



n 



5. Conclusion 



In this paper, we have shown that there exists a probabilistic counterpart to the classical operational calculus. 
A further direction of research is the extension of this approach to the multivariate case. 
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